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1. INTRODUCTION 

The studies of strongly interacting fermion systems are of considerable interest 
in solid state physics, nuclear physics, astrophysics, and so on. It is well known that 
many of the basic properties result already from the notion of the quasiparticles, 
at temperature T = 0, characterized by the step-like occupation function n(p, a) = 
np(p, cr) = 6(pf — p)i so that the number density, p = J2 a I n (P' cr )^P/(2 7r ) 3 , deter- 
mines the Fermi momentum p^ F = 37r 2 p [1]. 

In real physical systems, presence of various interactions between particles (by 
Coulomb forces, mesonic forces, van der Vaals forces, etc.) can essentially change the 
ground state energy and can rarely diversify the phenomenology of a Fermi system, 
preserving it as a Landau Fermi liquid. The Landau theory of normal Fermi liquids 
is based on the notions of the quasiparticles and the amplitudes F which characterize 
the effective interaction between the quasiparticles [1]. It permits to describe all the 
observable physical properties in terms of the first and second variational derivatives 
of the energy functional E[n(p, a)] with respect to the occupation numbers n(p, a): 
e(p, a) = 8E/8n{p, a) and F(p, a, p', a') = 5 2 E/5n(p, a)5n(p', a'). As a result, the 
Landau theory has removed high energy degrees of freedom and kept a sufficiently 
large number of relevant low energy degrees of freedom to treat liquid's low energy 
properties. Usually, it is assumed that the breakdown of the Landau theory is defined 
by the Pomeranchuk stability conditions and occurs when the Landau amplitudes 
being negative reach its critical value [2]. The new phase at which the stability 
conditions are restored can in principle be again described in the framework of the 
Landau theory. But it is not a general case. The most evident example is the 



superconducting transition under arbitrarily weak attraction between fermions with 
opposite spins a and opposite momenta close to the Fermi momentum p-p [3] . But one 
can also search for alternative scenarios of such restructuring produced by repulsive 
interactions, if strong enough. The quasiparticle dispersion law 

( \ 6E m 

e(P><r) = T-7 r 1 

Sn{p,a) 

satisfies the stability criterion: vf = [de(p, cr)/dp] p=PF > 0. Then the minimum of 
the energy, 

6 ( E -^ I n(p ' a) (^) = / (e(P^)-/*)MP^)(^a=0, (2) 

is reached (although with the chemical potential \i generally different from p F /2m) 
still for n = uf- This quasiparticle occupation should not be confused with the 
renormalization of real particle occupation function [4]. 

However one can also search for more general solutions of Eq. (2), n ^ uf- 
Regarding the Pauli principle restriction < n(p, a) < 1, these general solutions 
must satisfy the equation 

n (P, o-) [1 - n (P, o-)] [e(P» o-) - lA = °- ( 3 ) 

The possibility for the last factor in the l.h.s. of Eq. (3) to be zero implies a 
possibility for the so-called Fermi condensate (FC), manifested by a flat dispersion 
law: e(p, a) = fj,, and a fractional occupation < n(p, a) < 1 within a certain finite 
spherical layer pi < p < pf in the momentum space. Otherwise, the occupation 
function can only take the values or 1, so that e(p, a) < jU, n(p, a) = 1 at p < pi 
and e(p, a) > [i, n(p, a) = at p > p/ [3,5-7]. But the full set of alternatives for the 
ground state includes, besides FC and Fermi ground state (FGS), also the so-called 
multiconnected Fermi states (MFS). This possibility was first recognized yet in the 
Hartree-Fock framework [8,9]. We will study such a possibility using model forms for 
the interaction: F(p, a, p', a') = g a ,a'U aia '{q), with only relevant dependence on the 
momentum transfer q = |p — p'| and some coupling constants g a ,a'- 

Let us restrict consideration to uniform, isotropic, and spin- symmetric Fermi 
systems, where e(p, a) = e(p) (this is also the scope in which all the known FC models 
were considered). Then any MFS presents a sequence of fully occupied areas (we call 
them icebergs) where e(p) < /U, separated by empty spacers where e(p) > fx. Hence it 
can be uniquely labeled by the number iV of interfaces between icebergs and spacers: 
MFS(iV). When this number is odd, the MFS(iV = 2n — 1) consists of n icebergs: 
P2i-i < p < P2i (i = 1 , ■ ■ ■ n), and n spacers: p2i-2 < p < P2%-\ (including a central 
spherical void < p < pi). For even numbers, the MFS(iV = 2n) contains n icebergs 
(including a central spherical island < p < pi) and n — 1 spacers. Interfaces pi are 
fixed by the equilibrium conditions 



s(pi) = Li, i=l,...N, 



(4a) 



and icebergs fit to the normalization condition: 

n 

E^-^-i) = 37r 2 p (46) 

i=i 

(for odd N, the analogue for even iV is evident). In this context, FGS can be identified 
as MFS(O). 

At least, the choice between the two different types of FGS restructuring with 
growing repulsion constant g is uniquely defined by the analytic properties of the 
model potential U(q), so that FGS— >FC scenario is only possible if this function (or 
its derivatives) have singular points at real axis in the complex (/-plane, otherwise 
the infinite sequence of transitions FGS^MFS(l) ^MFS(2) -> . . . takes place [10]. 

Despite such strict dichotomy between the two evolution routes of model fermion 
systems, several analytic and numerical studies [10,11] showed that their expected 
physical properties should become very close in this course. If the distinction between 
analytic and non-analytic potentials is not fully respected, one can even arrive at a 
wrong conclusion on possible MFS(iV) ^FC transition at some finite N [12]. In fact, 
as will be shown below, the FC state is only reached in the limit of MFS(iV — > oo) (for 
a given analytic U(q)). But aside of these formal issues, there remains a clear interest 
in practical conditions when the difference between the MFS and FC states cannot 
be ignored, permitting to use the MFS picture for explanation of various actual 
physical effects related to the behavior of heavy-fermion systems in magnetic fields 
[13,14]. The present work is just focused on a detailed analysis of common features 
between FC and advanced MFS(iV ^> 1) states. To this end, we use the analytic 
model U(q) = 1/ \/ q 2 + q$ which tends to the non-analytic FC model 1/q in the limit 
qo — > 0, to develop an effective approximation for the average occupation function 
in this limit and to establish the criteria on model parameters (g, qo) and external 
parameters (as temperature) for the two states to be practically indistinguishable. 



2. ANALYTIC SOLUTION FOR THE l/y/q 2 + q% MODEL 
Let us start from the model interaction 

F(p, g ,pV)= 9 2 . (5) 

v(p - P ) 2 + % 

This form implies spherically symmetric and spin independent occupation: n(p, a) — > 
n(p), and we can present the energy functional as a ID integral: 

E[n{p)} = J p 2 n{p) [eoip) + V{p)/2] dp. (6) 

The potential energy V(p) also enters the quasiparticle dispersion law: 



s(p,a) =£ Q {p) +V(p), 



(7) 



here So(p) = p 2 /(2m) is the kinetic energy of system's particle with the bare mass 
m. V(p) presents a function of p and a functional of n(p') given by a ID integral: 



V(p) = ^ j p' 2 n{p')v{p,p'W : (8) 

where the kernel v(p,p') results from integration of the model potential U(\p — p'|) 
in the solid angle between the vectors p and p': 

i rP+p' 

v(p,p') = ^ tU(t)dt. (9) 

2 PP J\p-p'\ 

In the limit qo — > 0, the model, Eq. (5), passes into the FC model £/fc(<?) = l/<? 
[5,7]. Recall briefly its main properties in the present notation. The kernel, Eq. (9), 
is here simply 

vfc(p,p') = l/max(p,p') , (10) 
At small enough g, when FGS still holds, this gives the explicit potential energy 

f 3^|/(2m)-^ 2 /(2m), p < p F , 

Vfc(p) = I (11) 
{ t,pl/(rnp) , p > p F . 

When £ = gm/(3n 2 ) — > 1, the dispersion in Eq. (7) vanishes for all occupied states, 
manifesting the onset of FC, which leads to the chemical potential //fc that is lower 
then the corresponding to FGS. For £ > 1, the FC occupation function is: n-Fcip) = 
^~ 1 0(pf—p) withp/ = Pf^ 1 ^ 3 , thus keeping the quasiparticle energy constant: e{p) = 
3py/(2m), within the whole occupied band p < Pf (pi = in the 1/q model). 

In the actual model, Eq. (5), the kernel takes a more complicate form 



v(p,p) = 



\/(p + p') 2 + q 2 ~ V(P-P') 2 + 



2 



2pp' 



and using it in Eq. (8) at small enough g when FGS holds, we calculate the potential 
energy produced by the fully occupied Fermi sphere of radius pp: 

V(p) = V PF (p) + V PF (-p), (12) 

where 



VWTWT7 f p2 + ql + p p-p\ _ M arcs . nh p + P 



V V 2 y 2 q 




Figure 1. Instabilities of Fermi states in the £/ a/c/ 2 + c/ 2 model at qo = O.lpp 
with growing £. a) Dispersion law e(0) tends to [i from below at £ — > £1 ~ 1.04293. 
b) A void is growing in the center of filled Fermi sphere at £1 < £ < £2, c) an island 
will appear in the center of the void at £ = £2 ~ 1.04356. 



Note that, unlike Eq. (11), the form of Eq. (12) defines analytic behavior of 
e{p) for all real p. In presence of finite "screening parameter" qo, the repulsion 
potential U(q) is reduced compared to Ufc(q), hence the critical value of £, when 
FGS fails, turns higher than 1. It is easy to see that V(p) is an even function of 
p with maximum at p = and satisfies the inequality: V(p) > V(0) + p 2 V"(0)/2. 
Hence the instability of FGS in this model occurs just in the center of the Fermi 
sphere, when e(0) — ► s(pf) = A* (Fig- l a )- As discussed below, this instability is the 
first in an infinite series, and the above condition involves the corresponding values 
of potential energy: 



V(0) = ^ (pF\Jpl + ql - c^arcsinh^ 



and 

V( PF ) 
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\pl + gl) ^ + q « - ^arcsinh^ F ( 90 ' ' 



Pf 2 qo \p F 

Using them in Eq. (1) we obtain the first critical value of £ as a function of qo: 



6(90) =[3— g 

IK-"" h i) + (s)' rl (i3) 

^0), 1+ (^) 2 lu(^) 

in the important limit qo <C p-p and as £((?o) ~ (Qo/pf) 3 at c/o ^> Pf- At £ — ■> £1(170)5 
the dispersion law in the center is characterized by a negative effective mass: m* = 



It behaves as 



l/e"(0) = m/[l — £i(<Zo)] (insert to Fig. la). When £ exceeds £i((/o) 5 a spherical 
void of small radius pi opens in the center, followed by the first iceberg pi < p < P2 
(Fig. lb). This corresponds to the FGS — > MFS(l) transition, and the values of p\ 
and P2 result from the equilibrium and normalization conditions, Eqs. (4). Their 
numerical solution shows that the void radius extends with £, and in this course the 
negative effective mass diminishes and changes to positive. Eventually this again 
leads to an instability in the center: e(0) — > e(pi) = e(j>2) = A*, when £ attains the 
next critical value £2(90) > £i(<?o) (Fig- 1°)- But in this case a filled spherical island 
(second iceberg) emerges in the center of the void, and the system passes to MFS(2). 
With growing £, the alternating formation of voids and islands in the center continues 
infinitely, and the number of icebergs increases very rapidly, making the numerical 
analysis of exact conditions, Eqs. (4), practically impossible. 

However, an effective asymptotical treatment can be proposed for advanced 
MFS(iV 3> 1), suggested by the general £-g phase diagram. Considering the mono- 
tonically growing function £i(qo), Eq. (13), and the fact that £jv+i(<7o) > £n(Qo), 
we reasonably expect that MFS(iV ^> 1) are reached with decreasing screening pa- 
rameter q at fixed £ > 1. The averaged characteristics of such state should be close 
to the FC characteristics for the same £. Below we develop an effective description 
of MFS(iV 3> 1) using this closeness and prove by the obtained results that such 
approximation is asymptotically correct. 

3. EFFECTIVE DESCRIPTION OF MFS(iV > 1) 

We expect that, for an MFS(iV 3> 1) at given £ > 1 and q <C p/, the average 
occupation function 



(at pf ^> A qo), is close to the corresponding FC function npcip) = £ _1 #(p/ ~ P)i 
constant over the whole range [0,p/]. Therefore we model the "microscopic" MFS(A^) 
function n(p) by a certain effective function n e s(p). In vicinity of any < p < pf 
it consists in a sequence of icebergs of width as 5(p)/£ with spacers of width m 
(1 — l/£)5(p) (Fig. 2). The local period 8{p) is supposed to be a slow function of p 
(such that S'(p) 1) chosen in order to satisfy the equilibrium conditions, Eqs. (4), 
for the effective dispersion law e e s (p) . Then the overall multiplicity is asymptotically 
given by 




(15) 




(16) 



We present the effective potential energy as: 



N 




(17) 



where ith iceberg with the center in Pj (Fig. 2) contributes by: 




(18) 



Si = 5 (Pi). The most important step consists in separating from V e g(p) the FC 
function Vfc(p), so that e e ff(p) = /Upc + SV(p). Then the small functional SV(p) = 
^eff(p) — Vfc{p) can be naturally approximated by the linear terms in (effectively) 
small differences 5n(p) = n e s(p) — npc(p) and Sv(p,p') = v(p,p') — v^cip, p')'- 



where 



SV(p) = SV 1 (p) + SV 2 (p), 



SV 1 (p) 
SV 2 (p) 



3 
m 

M 

m 



p' 5v(p,p')dp' , 

fP 2 f fPf 

I p Sn(p')dp' + / p'Sn(p')dp' 
o Jp 



(19) 



A direct calculation shows that 8V\(p) is a negative monotonously growing function 
such that 

and Vi(p) ^ at p ^> p/. When integrating in SV 2 (p) over each interval Si centered 
in Pi, we first suppose that the width of ith iceberg is exactly Under this 

condition, called the geometric equilibrium (Fig. 2a) , the simple rules hold: 



<5,/2 
Si/2 



/Oil A 
Sn(P l + x)dx = 0, 
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5n(Pi + x)xdx = , 



(20) 



Sn(Pi + x)x 2 dx = S] 



; 2~f 
2^ 



At a special choice £ = 2 1 / 3 , also the last line in Eq. (20) turns zero, and the 
contributions to V 2 (p) from all the intervals vanish, except for that containing p 
itself. Adopting the above choice, we obtain the only contribution to SV 2 (p) from the 
interval Pi — Si < p < Pi + Sf. 



V oac (p) 
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(Pi + x)dx . 



(21) 



The totality of such contributions defines a continuous positive function, reaching 
zero between the icebergs: V osc (Pi ± Si/ 2) = 0, and maximum at their centers: 



V oac (Pi) = \U(p) 



35 2 (p)(tj-l) 
4m£ 2 




Figure 2. a) Condition of geometric equilibrium between icebergs and spacers, b) 
Physical equilibrium between iceberg distribution and dispersion law at slow variation 
of period. 



so that its fast oscillating behavior can be modelled by a simple expression: 

2np 



V oac {p) = 



1 + cos 



8(p) 



(22) 



However, the equilibrium for MFS requires that the minima of £ e ff(p) be located 
within icebergs and its maxima within spacers, opposite to what happens with V osc (p) 
at geometric equilibrium (Fig. 2a). This misfit can be corrected with small correlated 
shifts of iceberg boundaries from the geometric equilibrium, described by a finite 
gradient 5'(p). Then non-zero contributions to 5V2(p) appear from all the intervals, 
giving rise to its slow component V s \ w(p) ~ S' (p)p 2 j- / m while the fast component 
K)sc(p) only obtains a certain phase shift (see below). For this situation, we present 
the chemical potential as 

H = li fc + SV t (p) + Kiow (p) + V b (p) , (23) 

where Vb(p) ~ V max (p) is the value of V osc (p) at iceberg boundaries. The constancy 
of fx is mainly controlled by the condition that V s \ w(p) scales with 5Vi(p), leading to 
the estimate 

y(p)~4 ln — • (24) 
P f Qo 

On the other hand, the extrema of V osc (p) estimated, e.g., from differentiation of Eq. 
(22), are shifted due to 8'{p) ^ by 

5\p)5'{p) 



5(p) -pS'(p) 



and they can reach physical equilibrium positions (Fig. 2b) if this shift is ~ 5(p)/2. 
This implies the condition on the gradient 5'(p) ~ S(p)/p, agreeing with the initial 



assumption S'(p) <C 1. Then the comparison with Eq. (24) gives an estimate for the 
slow function 

*(p)~p4ln^. (25) 

This approximately linear behavior is valid at p ^> qo and Pf — p ^> qo, otherwise it 
is limited by the edge effects, the limiting values being 5 min ~ (<?o/P/) m (P//?o) and 
^max ~ (Qo/Pf) m (P//"?o)- The oc 1/p decrease of the density of icebergs, defined by 
Eq. (25), reflects the outward pressure on them from the center. Then the estimate 
for total multiplicity follows from Eq. (16): 

W „ *L^'*! = 4- (26) 
Qo HPf/Qo) J qo P % 

Also, using Eq. (25) in Eq. (22), we confirm that V& is negligible in Eq. (23). 
Correspondingly, the chemical potential /i is shifted downwards with respect to the 
FC value /ipc by 

A^^ln^-^lniV, (27) 
m q N 

and, since the deviations ~ V^ sc <C A/j, the energy difference between FC and 
MFS (A 7 ") is orrectly estimated by Eq. (27). We notice that all the considerations 
after Eq. (23) are equally valid for the general situation, £ ^ 2 1 / 3 . 

The asymptotic behavior of MFS(iV) resulting from Eqs. (25-27) depends only 
weakly on the coupling constant £. However this almost universal regime is only 
reached for N above some crossover value N* which essentially depends on £. This 
relates to the fact that multiplicity A" for given £ grows with 1/qo, but this process 
starts from 1/qo = l/q% such that £ = £i (<?£), Eq. (13). Hence the universal regime 
is reached already at 1/qo > l/q% and corresponds to iV > iV* ~ (pf / q^) 2 . Thus, for 
i - 1 « 1, we estimate from Eq. (14) A^* ~ ln[l/(£ - !)]/(£ - 1). 



4. CONCLUDING REMARKS 

At T = 0, our treatment has shown that the energy difference between FC and 
MFS(iV) tends to zero when iV overcomes the crossover value N* , N > N* . For MFS, 
the effective mass M* is defined by the dispersion near the boundaries of icebergs: 

Hjf* P P mN iOQ\ 

M ~dI^ld-p~ m W)~^ (8) 

It is seen from Eq. (28) that at iV ^> A^* the effective mass can be extremely large, 
though finite even at T = 0. Thus MFS represents a heavy-fermion system which can 
be treated as Landau Fermi liquid. At finite temperatures, the system persists to be a 
Landau Fermi liquid, but there is a crossover temperature T* at which the difference 
between FC and MFS vanishes. To calculate T*, we observe that the single-particle 
spectrum should not be altered when the temperature reaches T* . The effective mass 
of a system with FC at finite T is given by [15] 



(29) 



Upon comparing Eqs. (28) and (29), we obtain 



Pf-Pi InN 

T ~ Tf ^ — w (30) 

where Fermi temperature Tp = p\l2k^,m. At T > T*, the system comes into the 
state with the effective mass, Eq. (29), and the difference between the MFS and FC 
states is eliminated. Hence, Eq. (30) defines the crossover temperature T* at which 
the system is caused to pass from a Landau Fermi liquid to a strongly correlated 
Fermi liquid with FC and temperature dependent effective mass. A more detailed 
analysis of the physical properties of this and other model Fermi systems and their 
possible relation to the observed effects in real systems will be the topic of forthcoming 
studies. 
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